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Two fossil groups of galaxies at z ~ 0.4 in the COSMOS: 
accelerated stellar-mass build-up, different progenitors 
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ABSTRACT 

We report on two fossil groups of galaxies at redshifts z = 0.425 and 0.372 discovered 
in the Cosmic Evolution Survey (COSMOS) area. Selected as X-ray extended sources, 
they have total masses (M 200 ) equal to 1.9 (±0.41) x 10 13 and 9.5 (±0.42) x 10 13 M , 
respectively, as obtained from a recent X-ray luminosity-mass scaling relation. The 
lower mass system appears isolated, whereas the other sits in a well-known large-scale 
structure (LSS) populated by 27 other X-ray emitting groups. The identification as 
fossil is based on the z-band photometry of all the galaxies with a photometric redshift 
consistent with that of the group at the la confidence level and within a projected 
group-centric distance equal to 0.5i?2oo, an d «ab ^ 22.5 mag limited spectroscopy. 
Both fossil groups exhibit high stellar-to-total mass ratios compared to all the X-ray 
selected groups of similar mass at 0.3 ^J z ^ 0.5 in the COSMOS. At variance with the 
composite galaxy stellar mass functions (GSMFs) of similarly massive systems, both 
fossil group GSMFs are dominated by passively evolving galaxies down to M stars ~ 
10 10 Mq (according to the galaxy broad-band spectral energy distributions). The 
relative lack of star-forming galaxies with 10 10 $J M 8tars ^ 10 11 Mq is confirmed by 
the galaxy distribution in the b — r vs i colour-magnitude diagram. Hence, the two 
fossil groups appear as more mature than the coeval, similarly massive groups. Their 
overall star formation activity ended rapidly after an accelerated build up of the total 
stellar mass; no significant infall of galaxies with M stars ^ 10 10 Mq took place in the 
last 3 to 6 Gyr. This similarity holds although the two fossil groups are embedded in 
two very different density environments of the LSS, which suggests that their galaxy 
populations were shaped by processes that do not depend on the LSS. However, their 
progenitors may do so. We discuss why the late merging of a compact group is favoured 
over the early assembly as a formation scenario for the isolated, low-mass fossil group. 

Key words: X-rays: galaxies: clusters - galaxies: groups: general - galaxies: evolution 
- galaxies: luminosity function, mass function. 



1 INTRODUCTION 

In early numerical simulations, the merging of a compact 
Currently guest astronomer at the MPE galaxy group was found to lead to the formation of a single 
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elliptical galaxy (hence a "fossil group") in a few billion 
years (Barnes 1989JJ. The final product retains an X-ray 
emitting halo of hot gas (Ponman & Bertram 1993). The 
whole scenario was proved to be plausible by the discovery 
of the fossil group RX J1340.6+4018 (Ponman et al. 1994). 

The standard observational definition of a fossil system 
is based on the existence of extended X-ray emission with 
a bolometric luminosity Lx > 10 42 h^g erg s _1 and of an 
R-band magnitude gap Am^ > 2 between the brightest 
and second brightest members within 0.5-R20CU (Jones, Pon- 
man & Forbes 2000). In theoretical studies, the analogous 
definition of a fossil system is less straightforward (see e.g., 
D'Onghia et al. 2005 for a discussion). 

On the basis of the existing observations, mostly lim- 
ited to the local Universe, fossil groups are expected to be 
a considerable population among the gravitationally bound 
systems of galaxies (Vikhlinin et al. 1999; Romer et al. 2000; 
Jones et al. 2003; La Barbera et al. 2009; Voevodkin et al. 
2010). For instance, their estimated spatial density is equal 
to 2.83 x 10 -6 /175 Mpc -3 for an X-ray bolometric luminosity 
above 0.89 x 10 42 /if 5 2 erg s" 1 (La Barbera et al. 2009). X- 
ray properties and dark- matter (DM) content of fossil groups 
are comparable to those of bright groups and poor clusters 
of galaxies, with total masses of ~ 10 13 -10 14 Mq, whereas 
their brightest members are typically as luminous as a clus- 
ter cD galaxy (M R < -22.5 + 5 log h 50 ). 

In recent high-resolution hydrodynamical cosmological 
simulations in a A cold dark matter (ACDM) universe, fos- 
sil groups have already assembled half of their final mass at 
redshifts z Js 1, and subsequently they typically grow by mi- 
nor mergers only, whereas non-fossil systems of similar total 
mass on average form later (D'Onghia et al. 2005; Dariush et 
al. 2007). The early assembly leaves sufficient time for galax- 
ies with a luminosity L ^ L* (assuming a Schechter 1976 
luminosity function) to merge into the central galaxy by dy- 
namical friction, which yields the magnitude gap defining 
fossil groups. This "early formation" scenario is alternative 
to the compact group merging scenario (Barnes 1989). 

However, no difference seems to exist in the distribution 
of neighboring faint galaxy density excess, distance from the 
red sequencqj, structural parameters, or stellar population 
properties between the brightest central galaxies (BCGs) in 
25 fossil groups at z ^ 0.1 and a sample of bright field ellip- 
ticals with a similar distribution in optical luminosity and 
X-ray properties (La Barbera et al. 2009). These authors 
argue that this evidence undermines the idea that fossil sys- 
tems represent the "end product" of a physical mechanism 



1 These numerical simulations adopted a 1:4 ratio of luminous to 
dark matter and assumed that individual galaxies had dark halos. 

2 -Ra is the radius within which the total mass density of a 
group/cluster is equal to A times the critical density of the Uni- 
verse (p c ). Correspondingly, A/a = A p c (z) (4-7r/3) i? A is the 
total mass enclosed by Ra- 

3 In a given colour-magnitude diagram, the red sequence is the 
locus populated by galaxies dominated by old, passively evolving, 
stellar populations (see Arimoto & Yoshii 1987). This locus is 
particularly tight if the two broad-band filters that define the 
colour bracket the so-called 4000 A break. Among coeval galaxies, 
the old, passively evolving ones exhibit the reddest colours, if dust 
attenuation is neglected (see Pierini et al. 2004 for a discussion). 



that occurred at high redshift; it rather suggests that fossil 
groups are the final stage of mass assembly in the Universe. 

Whether either interpretation is true or fossil groups 
correspond to a transient phase in the life of a galaxy group 
(von Benda-Beckmann et al. 2008; Dariush et al. 2010; Cui 
et al. 2011), whether they form a distinct family of groups 
(D'Onghia & Lake 2004; Khosroshahi, Ponman & Jones 
2007) or not (Sales et al. 2007; Zibetti, Pierini & Pratt 2009; 
Sun et al. 2009; Democles et al. 2010; Voevodkin et al. 2010; 
Cui et al. 2011), for several reasons it is important to ex- 
tend the present knowledge of fossil systems beyond the local 
Universe and towards lower total masses. 

If undisturbed for a very long time, fossil groups may 
represent the ultimate example of hydrostatic equilibrium 
in viriliazed systems (Vikhlinin et al. 1999). This conclusion 
is supported by a recent detailed study of two such systems 
based on high-quality XMM-Newton spectro-imaging data 
(Democles et al. 2010). 

In addition, fossil groups could represent the least mas- 
sive, baryonically closed systems (Mathews et al. 2005) , their 
baryon mass fraction being equal to the cosmic mean valuqj. 
If so, they could form a reference sample with a recent, quies- 
cent merging activity with which one can investigate the ori- 
gin of the average baryon deficiency within -R500 (~ 0.77?20o) 
found in X-ray emitting groups at Z ^ 1 (Giodini et al. 2009, 
hereafter referred to as G09; see also Teyssier et al. 2011). In 
particular, two pair-matched samples of fossil groups with a 
range in total mass and redshift, one with radio-loud BCGs 
and the other with radio-quiet BCGs, could be used to esti- 
mate the "instantaneous" boost of the entropy in the central 
regions of the ICM that is due to radio-galaxy feedback (cf. 
Giodini et al. 2010). 

From simple considerations on the hierachical nature of 
structure growth via merging (Milosavljevic et al. 2006), fos- 
sil systems are expected to be more frequent (by a factor of 
ten) at the group mass-scale than at the cluster mass-scale 
(~ 10 14 -10 15 M ). This is confirmed by the most recent 
study of fossil groups based on hydrodynamical resimula- 
tions of five regions in the Millennium Simulation and on 
galaxy formation models (Cui et al. 2011 and references 
therein). Therefore fossil groups may also help in under- 
standing the properties of more massive bound systems at 
z = 0, especially if it turns out that they foster the forma- 
tion of the most massive progenitor of the BCG of a cluster, 
as suggested by J 0454—0309 at Z — 0.26 (Schirmer et al. 
2010; cf. Barnes 1989 and references therein). 

During our ongoing studies of the X-ray selected groups 
in the 2 deg 2 area of the Cosmic Evolution Survey (COS- 
MOS; Scoville et al. 2007), we have discovered two fossil 
groups at z = 0.372 and 2 = 0.425. As demonstrated by 
G09 and Giodini et al. (2010), the multi-wavelength data 
set of the COSMOS offers a thus far unique opportunity to 
characterize X-ray selected groups of galaxies up to z ~ 1. 
This enables us not only to fully describe the galaxy popula- 
tions of the two fossil groups, but also to characterize their 



4 The cosmic mean value of the baryon mass fraction can be 
obtained e.g. from the five years measurements of the Cosmic 
Microwave Background (CMB) performed with the Wilkinson 
Microwave Anisotropy Probe (WMAP5; Dunkley et al. 2009). 
According to this estimate, it amounts to 0.171 ± 0.009. 
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surrounding large-scale structures (LSSs), and to speculate 
on their progenitors. In Sect. 2 we describe selection and 
properties of these two groups. Their galaxy stellar mass 
functions and colour-magnitude diagrams are determined 
in Sect. 3. Discussion of results and inference on evolution- 
ary scenarios appear in Sect. 4; conclusions are drawn in 
Sect. 5. 

As in G09, we adopt a ACDM cosmological model 
(fi m = 0.258, n A = 0.742) with H = 72 km s _1 Mpc -1 , 
consistently with results from WMAP5 (Dunkley et al. 2009; 
Komatsu et al. 2009). Thus, a redshift of 0.372 (0.425) cor- 
responds to a lookback time of 3.981 Gyr (4.409 Gyr), a 
luminosity distance of 1961.1 Mpc (2298.2 Mpc), an angu- 
lar diameter distance of 1041.8 Mpc (1131.8 Mpc), and a 
scale of 5.051 kpc/" (5.487 kpc/"). 

Hereafter magnitudes are expressed in the AB system 
unless otherwise noted. 



2 SAMPLE SELECTION AND PROPERTIES 

From the composite mosaic of the XMM- Newton and Chan- 
dra X-ray data mapping the COSMOS 2 deg 2 area, extended 
sources (i.e., groups and clusters) were detected down to 
a flux limit of 10 -15 erg s -1 cm -2 (Finoguenov et al. in 
preparation; see also Finoguenov et al. 2007). G09 selected 
91 groups/poor clusters at 0.1 < z < 1, which have an 
X-ray flux estimate at a significance better than 3a, after 
point-like source subtraction (see Finoguenov et al. 2007). 
For these groups, the membership of a galaxy (through a 
precise photometric redshifflj) and the total stellar mass in 
galaxies together with its uncertainty were robustly deter- 
mined as described in G09. The total mass of each group was 
determined from the I/X-M200 relation established through 
weak lensing analysis by Leauthaud et al. (2010). Uncertain- 
ties on M200 are statistical only (i.e., they arise from the 
L1X-M200 relation). The effects of systematics are estimated 
to be lower than the statistical uncertaint}|j. 

We searched for fossil groups in the G09 sample. As in 
Giodini et al. (in preparation), we defined as member can- 
didates of a group those galaxies within a projected group- 
centric distance of -R200, and with a photometric redshift 
consistent with the photometric redshift of the group within 
two times the uncertainty on the photometric redshift of the 
group. Then we inspected the brightest member-candidate 
galaxies within a projected group-centric distance equal to 
0.7i?2oo to minimize the potential bias produced by uncer- 
tainties on i?200- As for ranking in luminosity, the avail- 
able i-band photometry was considered, since the reference 



5 For galaxies with i ^ 25 mag at the redshift of the two fossil 
groups under study (i.e., z m 0.4), the l<r error on the photometric 
redshift is equal to 0.03 X (1 + z) = 0.042. 

6 Systematic effects tested were: errors on photometric rcdshifts, 
mis-centering, uncertainties in the mass— concentration relation. 
These amount to a few percents (see fig. 7 in Leauthaud et al. 
2010), while the statistical error is much larger. Removal of point- 
like sources can lead to underestimate the X-ray flux, since the 
scale at which they are removed matches with the scale of cool- 
cores. The correction, as detailed in sect. 4.2 of Leauthaud et al. 
(2010), amounts to 3—17% of the total X-ray flux. This correction 
was taken into account when establishing the L-^—M^oo relation. 



COSMOS photometric catalogue is extracted from imag- 
ing in this band (Ilbert et al. 2009). This choice avoids the 
use of template- and (photometric) redshift-dependent k- 
corrections as for the choice of considering rest-frame ri- 
band magnitudes. Furthermore, up to 2 = 1, the i-filter 
(observed frame) maps emission at rest-frame wavelengths 
redder than the 4000 A-break, which limits the impact of 
young stellar populations and, potentially, attenuation by a 
dusty interstellar medium. 

A number of groups were initially selected as poten- 
tial fossil systems, the magnitude difference between the 
BCG and any other bright member-candidate galaxy within 
0.7i?2oo being conservatively chosen to be larger than 1.5 i- 
mag. We set a larger search radius and a smaller magnitude 
gap than those required by the definition of fossil groups (cf. 
Sect. 1) in order to maximize the number of candidates in 
this first round of the search for fossil groups. Two groups 
were successively established as fossil thanks to the exist- 
ing spectroscopic information (Lilly et al. 2007), limited to 
sources brighter than i = 22.5 mag, which enabled a finer 
screening against projected foreground/background bright 
galaxies. 

On one hand, we flagged out galaxies classified as mem- 
ber candidates if they were located within -R200 from the 
X-ray centroid of either fossil group but exhibited a line-of- 
sight velocity in excess of 3000 km s -1 plus its 3er uncer- 
tainty with respect to the spectroscopically determined re- 
cession velocity of either group. The redshift of each group 
was computed using the median spectroscopic redshifts of 
the galaxies within -R200 through iterative removal of out- 
liers with a recession velocity beyond 1500 km s -1 with re- 
spect to the recession velocity of the group. The starting 
redshift was defined by the photometric redshifts of galaxies 
on the red sequence. Contamination was checked by run- 
ning a red-sequence-finding algorithm between redshifts 
and 2 and by using spectroscopic group catalogs (Knobel et 
al. 2009) . On the other hand, the available spectroscopy en- 
abled us to check, in particular, that bright galaxies within 
0.5i?2oo from the X-ray centroid of either fossil group were 
actually not dropped as member candidates because of an 
inaccurate photometric redshift. Additional considerations 
on size and luminosity were used to check the goodness of 
the photometric redshifts of relatively bright galaxies with- 
out spectroscopy (see CXGG 095951+0212.6 in Fig. 2). 

Figures 1 and 2 show cut-outs of the central regions 
of the two COSMOS fossil groups under study, which en- 
compass the area within 0.5i?2oo for each system. They re- 
produce X-ray surface brightness contours overlaid on the 
i-band HST/ACS images (Koekemoer et al. 2007) of the two 
galaxy systems. In particular, member candidates (i.e., pho- 
tometric members) are indicated together with all galaxies 
brighter than i = 22.5 mag having a spectroscopic redshift. 
The BCG of either group can be easily identified at the cen- 
ter of the X-ray emission region. 

One system, CXGG 095951+0140.8, has Am il2 = 
2.19 ± 0.014 mag (i.e., Ami 2 = 2.10 ± 0.02 mag as obtained 
from the difference between the magnitudes in the rest-frame 
r-band), R200 = 832 ± 19.6 kpc, M200 = 9.5 (±0.42) x 



7 The uncertainty of the high-quality (flags 3 and 4) zCOSMOS 
spectroscopic redshifts amounts to 110 km s — . 
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Figure 1. i-band image of CXGG 095951+0140.8, enclosing a circular area of radius 0.5i?200 centered on the X-ray position of the group 
(marked in red). Photometric members are marked with squares; green circles mark the spectroscopic members. Available spectroscopic 
redshifts of all objects in the region are given (yellow diamonds). Dashed contours indicate values of the X-ray surface brightness at the 



4<r, 7ct and 10<r levels after removal of point-like sources, where lc corresponds to 4.6 X 10 
band. 
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0.372, and six spectroscopic members. It 



is part of a large-scale structure (see fig. 1 in G09 and 
fig. 3 in Scoville et al. 2007) populated by 28 X-ray emit- 
ting groups distributed across the entire 2 deg 2 area of the 
COSMOS (corresponding to a cross size of about 25.5 Mpc 
at Z — 0.37). Furthermore, its dominant galaxy (with a 
rest-frame absolute magnitude Mi = —24.87) hosts a point- 
like radio source (see Giodini et al. 2010). The other fossil 
group, CXGG 095951+0212.6, has Am il2 = 2.35+0.014 mag 
(i.e., Anii2 = 2.32 ± 0.02 mag in the rest-frame r-band), 
#200 = 478 ± 54.4 kpc, M 200 = 1.9 (±0.41) x 10 13 M , 
z = 0.425, and eight spectroscopic members. It is isolated 
and its BCG has Mi = -23.87. 

Basic properties of the two fossil groups under study 
are listed in Table 1. In addition, we note that they pop- 
ulate the upper half of the distribution of the COSMOS 
X-ray selected groups in the (galaxy) stellar mass fraction- 
group total- mass diagram (G09, their fig. 5), where quan- 



tities are estimated at 0.7i?2oo- This is particularly true 
for CXGG 095951+0212.6. As from Gofl the stellar mass 
fraction is equal to 0.039 for the more massive fossil group 
(CXGG 095951+0140.8) and 0.143 for the less massive one 
(CXGG 095951+0212.6). Thus, the fraction of the stellar 



8 Giodini et al. (in preparation) use a new catalogue of galaxy 
stellar masses by O. Ilbert instead of the catalogue by Ilbert et al. 
(2009) used in G09. The difference between the catalogues rests 
on the way stellar masses are computed. In the new one, stellar 
mass estimates are based on fitting the spectral energy distribu- 
tion (SED) of a galaxy and not on the rest-frame K-band lumi- 
nosity of the latter; furthermore, the stellar initial mass function 
of Chabrier (2003) is used instead of a Salpeter (1955) one. For 
the purposes of our study it is important that no systematics dif- 
ferently affects such estimates in fossil and non-fossil groups, not 
that the absolute scale of the estimated stellar mass of a galaxy 
is correct. 
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Figure 2. The same as in Fig. 1 but for CXGG 095951+0212.6. 
3.5 X 10 -16 erg s _1 cm -2 arcmin -2 in the 0.5—2 keV band. 



Here the la value of the X-ray surface brightness corresponds to 



mass in galaxies decreases with increasing total mass for 
the two fossil groups as for all other coeval groups (G09). 

Furthermore, the more (less) massive fossil group ex- 
hibits a value of the total i-band-to-X-ray luminosity ra- 
tio at -R200 (Li/Lx) equal to 138 (305), after background 
subtraction including spectroscopic information (from Gio- 
dini et al. in preparation). Therefore CXGG 095951+0140.8 
(CXGG 095951+0212.6) exhibits a low (normal) value for 
its total mass with respect to the distribution of Li/Lx as a 
function of A/200 for the 51 groups at Z — 0.3-0.5 considered 
by Giodini et alQ 

With a total mass of 1.9 (±0.41) x 10 13 M , 
CXGG 095951+0212.6 is among the least massive fossil sys- 
tems selected in X-rays, whereas CXGG 095951+0212.6 is at 
the boundary between the group and cluster mass regimes 
(cf. Vikhlinin et al. 1999; Sun et al. 2009; Zibetti et al. 2009; 
Democles et al. 2010; Voevodkin et al. 2010). As a refer- 
ence, RXJ105453. 3+552102 at z = 0.47, with a total mass 



9 These COSMOS X-ray selected groups define a mass-complete 
sample, probing different large-scale density environments. 



of ~ 10 Mq (Aguerri et al. 2011), is among the most mas- 
sive fossil systems 10 !. 



10 RXJ105453. 3+552102 is not considered as a fossil system 
by Aguerri et al. (2011), since its magnitude gap is equal to 
1.92 ± 0.09 mag in the observed-frame r-band. The bright end of 
the galaxy luminosity function of this cluster is populated by old, 
passively evolving objects (Aguerri et al. 2011), so that the intro- 
duction of a fc-correction is expected not to significantly modify 
the magnitude gap in the rest-frame r-band. However, the fos- 
sil criteria fixed by Jones et al. (2003) are somewhat arbitrary 
and depend on the way magnitudes are measured (see Aguerri 
et al. 2011 for a discussion). In addition, the likely contribution 
from diffuse light to the BCG emission is hard to be fully ac- 
counted for (see Zibetti et al. 2009 for a discussion). Thus, we 
consider RXJ105453. 3+552102 as a fossil system as we did for 
RXCJ2315. 7-0222 at z = 0.026 with Am 12 = 1.88 ± 0.03 mag 
in the rest-frame Rc-band (Democles et al. 2010). 
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Table 1. Basic properties of the sample 



Denomination 


z 


M 20 o 
[10 13 M ] 


R200 
[kpc] 


Am il2 
[mag] 


Ami2 
[mag] 


N 


Nphot 


^spcc 


CXGG 095951+0140.8 
CXGG 095951+0212.6 


0.372 
0.425 


9.5 ± 0.42 
1.9 ± 0.41 


832 ± 19.6 
478 ± 54.4 


2.19 ± 0.014 
2.35 ± 0.014 


2.10 ± 0.02 
2.32 ± 0.02 


261 
68 


140 

29 


6 

8 



Denomination, spectroscopic redshift (z), total mass (M200); size (-R20o)i observed i-band and rest-frame r-band 
magnitude gaps between the first two brightest member galaxies (Am;j 2 and Ami2, respectively), total number of 
galaxies brighter than the completeness magnitude (i = 25.37 mag and i = 25.17 mag for CXGG 095951+0140.8 
and CXGG 095951+0212.6, respectively) within -R200 (^V)i numbers of expected (i.e., the photometric members 
after statistical subtraction of field galaxies) and spectroscopic member galaxies within -R200 (^Vphot an d N spec , 
respectively) are listed for each fossil group. 



3 RESULTS 

3.1 Galaxy stellar mass function 

We compute the galaxy stellar mass function (GSMF) of the 
two fossil groups in bins of 0.5 dex in stellar mass (M stars ). 
Only member-candidate galaxies with a stellar mass above 
the completeness mass at the redshift of each group are con- 
sidered. The completeness mass is equal to 8.3 x 10 8 Mq or 
1.25 x 10 9 Mq for, respectively, all member-candidate galax- 
ies of CXGG 095951+0212.6 or those among them that are 
classified as passively evolvin g 11 ! . It is equal to 3 x 10 8 Mq 
or 5 x 10 s Mq for, respectively, all member-candidate galax- 
ies of CXGG 095951+0140.8 or the passively evolving ones 
among them. 

We applied a correction for the expected number of 
interlopers, based on the statistics of the galaxies outside 
X-ray emitting bound structures in the same photometric- 
redshift slice, which define the coeval field of either group. 
This holds whether or not galaxies are additionally selected 
according to their star formation activity. In each case, the 
field GSMF was normalized by the group-to-survey volume 
ratio and then subtracted to the GSMF of the member- 
candidate galaxies to obtain the GSMF of the (bona fide) 
member galaxies of each fossil group. Counts per stellar mass 
bin are normalized to the total mass M200 of each fossil 
group, to enable comparison with the composite GSMF of 
the 40 (11) COSMOS X-ray selected groups at 0.3 < z sC 0.5 
with 2 x 10 13 < M200 < 5 x 10 13 M (M200 > 5 x 10 13 Mq) 
from Giodini et al. (in preparation). These two composite 
GSMFs are computed in an analogous way, but include the 
number counts from all selected groups with similar masses 
to that of either fossil group; they are normalized to the sum 
of the total masses of the individual groups. 

Figure 3 shows the normalized GSMFs of the two fos- 
sil groups and the normalized composite GSMFs of the two 
reference samples of X-ray selected groups at 0.3 ^ z ^ 0.5 
for all (bona fide) member galaxies and for those clas- 
sified as passively evolving. Comparison of these normal- 
ized GSMFs reveals that about 80% of the galaxies with 



Al 



> 10 1 



Mq are passively evolving in the two fossil 
groups at z ~ 0.4. Conversely, passive galaxies start dom- 
inate the GSMF of the coeval, mass-complete sample of 



X-ray selected groups at much higher stellar masses, with 
a fraction of 80% passive galaxies being reached only at 
M star S ^ 10 n Mq Thig hoMs although the f rac tion of 

passively evolving galaxies increases towards larger stellar 
masses across the entire range in stellar mass. 

Interestingly, half of the galaxies with M stars ~ 
10 11 Mq do not show evidence of ongoing star- formation 
activity for similarly massive, optically selected groups at 
0.3 < z sC 0.48 (Wilman et al. 2008 

If the fraction of star-forming ("blue") galaxies is a 
proxy for the dynamical youth of a galaxy system at fixed to- 
tal mass (e.g., Tully & Shaya 1984; Balogh, Navarro, & Mor- 
ris 2000; Ellingson et al. 2001; Zapata et al. 2009; Braglia et 
al. 2009), the previous difference suggests that the two fossil 
groups assembled their DM haloes earlier than the average 
system among the similarly massive groups at 0.3 ^ z ^ 0.5. 



3.2 Colour— magnitude diagram 

The dominance of passively evolving or quiescent galaxies at 
M stars ^ 10 10 Mq is confirmed, in a more direct way, by the 
b—r vs i colour-magnitude diagram of either fossil group (see 
Fig. 4). This diagram is obtained using Subaru photometry 
(see Ilbert et al. 2009). For z ~ 0.4, the b — r colour nicely 
brackets the redshifted spectral feature at ~ 4000 A (rest 
frame), whose amplitude increases together with the time 
since the last main episode of star formation but also with 
increasing metallicity (e.g. Worthey 1994). 

Both fossil groups exhibit a well defined red-sequence, 
and a lack of galaxies within 0.5i?2oo that are bluer than this 
locus by ^ 0.5 mag. Most of these blue galaxies are at least 
~ 4 magnitudes fainter than the BCG of either fossil group 
in the i-band. Moving out to R200, the bona fide member 
galaxies with a b — r colour at least 0.5 mag bluer than the 
BCG one are still a few. This confirms that most of the 
luminous galaxies within the bound region of either fossil 
group are old and passively evolving or quiescent, modulo 
dust attenuation (cf. Wilman et al. 2008). 

The dominance of galaxies with a very low star- 
formation activity among those that are about two to four 
magnitudes fainter than the BCG is confirmed for the 
fossil groups RXCJ2315. 7-0222 and RXCJ0216. 7-4749 



11 COSMOS galaxies are classified as passively evolving or star 
forming on the basis of the template best-fitting their individual 
SEDs. In particular, the former galaxies have star formation rates 
lower than 10 -2 Mq yr _1 . 



12 In Wilman et al. (2008), passive galaxies are defined on the 
basis of their Spitzer-IRAC mid-infrared colours, stellar masses 
are derived from rest-frame K-band luminosities (in analogy with 
Ilbert et al. 2009). 
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Figure 3. The galaxy stellar mass functions (GSMFs) of the two fossil groups CXGG 095951+0212.6 (left) and CXGG 095951+0140.8 
(right). Each panel shows the GSMF for all the (bona fide) member galaxies (black circles) and the GSMF for the subsample of passively 
evolving (bona fide) member galaxies (red symbols), together with their statistical uncertainties. Conversely, the vertical black and red 
lines represent, respectively, the mass completeness limit for all the galaxies at the same redshift as the group and for the passively 
evolving ones among them. In each panel the black and red dashed curves represent, respectively, the GSMFs for all the galaxies in 
the coeval field and for the passively evolving ones among them. To obtain the GSMF of all the (bona fide) member galaxies of either 
group, the corresponding field GSMF was normalized by the group-to-survey volume ratio and then subtracted to the GSMF of all the 
candidate members. The GSMF of the passively evolving (bona fide) member galaxies of either group was obtained in an analogous way. 
Furthermore, each panel shows the composite GSMFs of the coeval, X-ray selected groups for all their (bona fide) member galaxies (yellow 
area) and the passively evolving ones among them (red area) together with their statistical uncertainties (Giodini et al. in preparation). 
Among the X-ray selected groups at 0.3 ^ z ^ 0.5 in COSMOS, there are 40 systems as massive as CXGG 095951+0212.6 (i.e., with 
2 X 10 13 < A/200 ^ 5 X 10 13 M Q ) and 11 systems as massive as CXGG 095951+0140.8 (i.e., with M 2 oo > 5 X 10 13 M ). 



(Democles et al. 2010), AWM4 (Zibetti et al. 2009; see fig. 2 
in Koranyi & Geller 2002), RX J1416.4+2315 (see fig. 2 
in Cypriano, Mendes de Oliveira, & Sodre 2006), and the 
fossil group candidate RX J1520.9+484cE] (Eigenthaler & 
Zeilinger 2009). However, in the literature there is no com- 
parative study of the galaxy populations of fossil and non- 
fossil systems of similar masses down to such fainter magni- 
tudes, based on multi-band photometry or spectroscopy. 

Finally, we note that two member-candidate galax- 
ies with i ~ 23.4 and i « 24.1 exhibit colors that are 
much redder than average color of the red-sequence of 
CXGG 095951+0140.8 (6 - r « 2.3 and b - r fa 2.9, re- 
spectively, vs b — r ~ 1.7). These two very faint and red 
galaxies (about 6.3-7 i-mag fainter and 0.4-1 mag redder 
in b — r than the BCG) are located at projected distances 
of 110 and 330 kpc from the BCG (or the cluster cen- 
ter), respectively. Member-candidate galaxies with analo- 



gous photometric properties exist in the fossil group can- 
didate RX J1548.9+085£3 (Eigenthaler & Zeilinger 2009). 



4 DISCUSSION 

4.1 Characteristic timescales 

In order to understand why the bulk of the galaxies with 



M 



> 10 1 



10 Mq are dominated by old, passively evolving 
stellar populations in the two fossil groups at z ~ 0.4, we 
need to estimate if there was enough time for the last formed 
stars to age once their host galaxy had reached the bound 
region of either group. Hence, we need to estimate the fol- 
lowing timescales: the time elapsed since half the DM halo of 
each fossil group was in place (i.e., the so-called "assembly 
time"), and the time elapsed since the last significant infall 
of blue galaxies with M stars ^ 10 10 Mq occurred. These two 



13 For this object the virial radius was estimated on the basis of 
a scaling relation (Evrard, Metzler & Navarro 1996). 



14 The fossil nature of this group comes under scrutiny however 
(P. Eigenthaler private communication). 
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Figure 4. b — r vs i colour-magnitude diagram (observed frame) for all member-candidate galaxies of either fossil group (black circles); 
galaxies within a group-centric distance of 0.5.R200 are marked by filled circles. Red open squares mark those photometric-member 
galaxies with a spectroscopic redshift inconsistent with that of either group (see Sect. 2). 



timescales need to be compared with the minimum life-time 
of the stellar populations characterizing old, passively evolv- 
ing (i.e., red-sequence) galaxies, which is about 1-2 Gyr. 

We now make the conservative assumption that the two 
fossil groups behave as other similarly massive groups ob- 
served at z ~ 0.4. Furthermore, we assume that groups do 
not significantly increase their total masses between z ~ 0.5 
and z w 0.4 15 l. in a statistical sense. By doing so, we can 
use results on the median formation redshifts of DM haloes 
obtained from the simulations of Gao et al. (2004) or dif- 
ferent analytical models that are reproduced in fig. 5 of 
Giocoli et al. (2007). Groups of similar total masses to 
CXGG 095951+0212.6 and CXGG 095951+0140.8 have as- 
sembled half of their DM halo at typical redshifts of ~ 1.2 
and ~ 1.05, respectively, which correspond to similar look- 
back times of w 4 Gyr at z — 0.425 and z — 0.372. 

According to the computations of dynamical friction 
and galaxy-merging timescales of Boylan-Kolchin, Ma & 
Quataert (2008, their fig. 1), galaxies with M stars ~ 10 10 
or 10 11 Mq need about 6 or 3 Gyr, respectively, to reach 
an orbit with maximum group-centric distance of 0.5i?200 
from 7?200 for the first time. Since the simulations of Boylan- 
Kolchin et al. mostly neglect the effect of the central galaxy 
in a given DM halo on the merging time-scale, and do not 
take into account the drag on ambient gas (e.g. Ostriker 
1999) , we infer that the epoch of last infall of blue galaxies 
with Af stars > 10 10 M Q can date back to less than 3-6 Gyr 
with respect to the lookback time of either fossil group. 

The time of last infall of blue, intermediate-mass galax- 
ies is comparable to the median time elapsed since the as- 
sembly of a group as massive as either fossil group, if not 
longer than that. Thus, we can assume that most of the 
galaxies with M stars 



^ 10 10 Mq have resided within the 



15 According to the adopted cosmological model (Sect. 1), this 
redshift range corresponds to a time interval of 0.7 Gyr. 



bound region of either fossil group since its assembly (cf. 
von Benda-Beckmann et al. 2008; Dariush et al. 2010). This 
"residence" time of m 4 Gyr is longer than the typical life- 
times of young stars, which are responsible for the bulk 
luminosity shortward of the 4000 A-break, when they are 
present. Hence, there is enough time to observe a domi- 
nating population of old, passively evolving galaxies with 
M stars v, 1Q io Mq in the twQ fossil groups at z ~ 0.4. This 

is even more the case if the assembly time of a fossil group 
is shorter than the median assembly time for all coeval, sim- 
ilarly massive groups (D'Onghia et al. 2005; Dariush et al. 
2007). Environmental effects could be responsible for the 
raise of the population of old, passively evolving galaxies in 
the two fossil groups (see Boselli & Gavazzi 2006 for a recent 
review) . 



4.2 Star formation history 

Member-candidate galaxies of either fossil group which are 
satellites at 2 w 0.4 could reach stellar masses of 1- 
8 x 10 10 Mq before their star formation turned to its end. 
We note that, in spite of the generally held, negative feed- 
back exerted by the group environment on the star formation 
activity of member galaxies, such an activity seems to be en- 
abled in 20%-50% of the galaxies with M stars - 10 11 M in 
groups at similar redshifts and with similar masses (Fig. 2; 
Giodini et al. in preparation; see also Wilman et al. 2008 
and references therein). 

It is easy to understand how tidal interactional (e.g., 
Spitzer & Baade 1951; Richstone 1976; Farouki & Shapiro 



16 Tidal interaction may have been active in a de-localized way, 
in particular during the process of rapid and violent formation of 
the dominant galaxy (D'Onghia et al. 2005; see also Pierini et al. 
2008 for an analogy), but, also, in a more localized way, since the 
early formation of the central dominant galaxy (see Elmegreen et 
al. 2000 for an analogy). 
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10 Mq. It is less certain how the same in- 
^j could have fostered the growth of more 



1981; Merritt 1983, 1984; Icke 1985; Miller 1986; Byrd & Val- 
tonen 1990; Valluri & Jog 1990) or harassment (Moore et al 
1996) could have inhibited the growth of satellite galaxies 
up to M stars ~ 10 
teraction processe; 
massive satellite galaxies (e.g., Keel et al. 1985; Kennicutt 
et al. 1987; Byrd & Valtonen 1990; Mihos et al. 1992; Hen- 
riksen & Byrd 1996; Lake, Katz & Moore 1998; Iono, Yun 
& Mihos 2004; Braglia, Pierini & Bohringer 2007). 

An earlier assembly time of fossil groups with respect 
to non-fossil groups of sinrilar masses (D'Onghia et al. 2005; 
Dariush et al. 2007) can justify the lower fraction of star- 
forming galaxies with 10 10 sC AP^™ ^ 10 11 M (see 
Sect. 4.1). As discussed by D'Onghia et al. (2005), this is also 
consistent with the larger X-ray-to-optical luminosity ratios 
of fossil groups with respect to similarly massive non-fossil 
groups in the local Universe claimed by Jones et al. (2003). 
In fact, in groups with earlier assembly times the hot gas 
was collected and started cooling earlier and at higher cen- 
tral densities. However, no difference in the scaling relation 
between X-ray and optical luminosities is reported for a sam- 
ple of seven fossil groups at z ^ 0.2 and total masses M200 
of a few times 10 14 Mq (Voevodkin et al. 2010), at vari- 
ance with previous results (Jones et al. 2003; Khosroshahi 
et al. 2007). As mentioned in Sect. 2, the X-ray-to-optical 
luminosity ratio is rather high for CXGG 095951+0140.8 but 
normal for CXGG 095951+0212.6, when compared to those 
of similarly massive groups at z — 0.3-0.5. 

Interestingly, there is growing evidence that an environ- 
ment characterized by high galaxy number density and low 
galaxy velocity dispersion, such as the one found in compact 
groups, plays a key role in accelerating galaxy evolution by 
enhancing star-formation processes in galaxies and favor- 
ing a fast transition to quiescence (Tzanavaris et al. 2010; 
Walker et al. 2010 and references therein). Such an environ- 
ment could represent the immediate progenitor of a fossil 
system. This is straightforward if such a progenitor is a com- 
pact group, and compact groups are interpreted as transient 
phases in the dynanrical evolution of looser groups (Barnes 
1989). Alternatively, an environment with the previous char- 
acteristics can be produced if the existence of fossil systenrs 
is primarily driven by the relatively early infall of massive 
satellites (D'Onghia et al. 2005; von Benda-Beckmann et al. 
2008). Tentatively this may be achieved through preferen- 
tial infall on orbits with typically lower angular momentum. 
Consistently, there is an indication of the presence of galax- 
ies in radial orbits in the external region of the fossil cluster 
RXJ105453.3+552102 (Aguerri et al. 2011). 

In a future study, we will explore if an efficient confine- 
ment of the hot, X-ray emitting gas (Mathews et al. 2005) 
could foster the rather large stellar mass fractions of the two 
fossil groups under study (cf. G09). 



4.3 Fossil group progenitors 

An earlier assembly time for fossil groups with respect to 
as massive non-fossil groups is consistent with the obser- 



17 Tidal interactions and harassment describe, respectively, pro- 
longed and collective impulsive interactions of type galaxy-galaxy 
or galaxy-gravitational potential well of the parent system. 



vational results discussed in Sect. 3; it is also predicted 
by simulations in a statistical sense (D'Onghia et al. 2005; 
Dariush et al. 2007). The evidence is a GSMF dominated 
by old, passively evolving systems at M stars ^ 10 10 Mq 
(Sect. 3.1) plus the corresponding existence of a tight red- 
sequence (Sect. 3.2). Since the two fossil groups at z w 0.4 sit 
in very different density environments of the LSS (Sect. 2), 
the conditions leading to the formation of a fossil system, 
in particular to the shaping of its galaxy population, do 
not appear to be intimately connected to the surrounding 
LSS. This is also consistent with simulations (von Benda- 
Beckmann et al. 2008; Cui et al. 2011). 

An early formation time was advocated by Aguerri et 
al. (2011) to explain the dynamics and luminosity function 
of the member galaxies as well as the photometric properties 
of the BCG for the relaxed, massive (M200 ~ 1 x 10 15 Mq) 
fossil system RX J105453.3+552102 at z = 0.47. On the 
other hand, a deep optical image of the BCG of the low 
X-ray luminosity fossil-group candidate RXJ1548. 9+0851 
(Santos, Mendes de Oliveira & Sodre 2007) clearly reveals 
the existence of many shells (Eigenthaler & Zeilinger 2009). 
According to Hernquist & Spergel (1992), many shell sys- 
tems among elliptical and spheroidal galaxies may have orig- 
inated through major mergers of disc galaxies with compara- 
ble masses rather than via minor mergers - between a mas- 
sive elliptical galaxy and a small disc companion (Quinn 
1984; Hernquist & Quinn 1988, 1989) or a small ellipti- 
cal/spiral conrpanion (Dupraz & Combes 1986) - or accre- 
tions (Thomson & Wright 1990; Thomson 1991). We specu- 
late that the relatively early infall of massive satellites can 
better justify the formation of a rather massive fossil group 
in a denser region of the LSS, like CXGG 095951+0140.8, 
whereas compact groups can more likely be the progeni- 
tors of observed low-mass and isolated fossil groups, like 
CXGG 095951+0212.6. 

On one hand, the LSS surrounding the progenitor of 
a fossil group could be responsible more for the growth of 
the total mass of the system than for the growth of the 
stellar nrass of its BCG (see Khosroshahi, Ponman & Jones 
2006). This is consistent with both the shallow dependence 
of the median luminosity of BCGs on the total mass of the 
parent halo observed in the local Universe (Yang, Mo & 
van den Bosch 2008) and the result from the hydrodynamic 
simulations of Cui et al. (2011) that X-ray bright groups are 
disproportionally rare in low-density regions. 

On the other hand, it is reasonable to assume that iso- 
lated fossil groups can be identified as such for longer times 
than non isolated ones of similar mass, because it is more 
difficult either that their magnitude gap is replenished by 
infalling galaxies, or that they fall into larger bound struc- 
tures. In a low density environment, a loose group has more 
chances to keep its identity and dynamically evolve into a 
fossil group after a compact group phase (cf. Barnes 1989). 
Either in the case of early assembly (D'Onghia et al. 2005; 
Dariush et al. 2007) or in the case of final evolution of a 
loose group (Barnes 1989), a rapid build-up of the stellar 
mass and the following quick end of the star-formation ac- 
tivity in most of the brightest member galaxies need to be 
invoked. Interestingly, this is what compact groups in the 
local Universe seem to have experienced (Tzanavaris et al. 
2010; Walker et al. 2010 and references therein). 

If two channels lead to the formation of a fossil group 



10 D. Pierini et al. 



and there is a sort of density environment bias on identi- 
fication and inference on the progenitor from observations, 
the stellar populations of the BCG of a fossil group do not 
necessarily have to be remarkably different from those of 
field, massive early-type galaxies of similar optical and X- 
ray luminosities at z ~ (cf. La Barbera et al. 2009). Fur- 
thermore, the size of environmental effects on typical red- 
sequence galaxies is still a controversial issue, and might be 
quite small at large stellar masses (cf. e.g., Thomas et al. 
2005, 2010; Pannella et al. 2009; Cooper et al. 2010). Hence, 
there might be less tension between our conclusions and the 
result of La Barbera et al. (20090. 

Future multi- wavelength studies of statistical samples of 
fossil groups with different masses and redshifts (e.g., Santos 
et al. 2007), located in different density regions of the LSS, 
will provide a test of our interpretation and hypotheses. 



4.4 Intracluster light 

Here we discuss if the contribution from low-surface bright- 
ness, diffuse stellar emission not associated with galaxies 
(the so-called intracluster light - ICL) could explain the 
rather low value of Li/Lx found in Sect. 2 for the more 
massive fossil group, CXGG 095951+0140.8. The impact on 
the total stellar mass budget of the ICL has already been 
pointed out in the literature, and with different scales, from 
observations (cf. e.g., Zibetti et al. 2005 and Gonzalez et al. 
2007). In simulations, most of the ICL seems to be associ- 
ated with the build-up of the BCG (e.g., Murante et al. 2007 
and references therein). The expected amount of ICL is un- 
certain as well as the observed one - 3% to 30% of all stars 
in groups at z = (e.g., Kapferer et al. 2010 and references 
therein) - and seems to be independent of the mass-scale 
of the galaxy system (e.g., Dolag, Murante & Borgani 2010 
and references therein). 

Groups of galaxies that formed early (as fossil groups, 
potentially) might be particularly effective at producing 
ICL, as interactions in the group environment begin to lib- 
erate tidal material. Galaxy systems as massive as the two 
fossil groups under study have smaller velocity dispersions 
than do massive clusters, meaning that interactions on the 
group scales can be slow and damaging as opposite to the 
impulsive high speed encounters, e.g. by harrassment, on 
the cluster scale. These slow tidal interactions can strip ma- 
terial from rotating disc galaxies into loosely bound tidal 
structures, forming long tails of stars and gas (D'Onghia et 
al. 2009). It is then interesting to estimate this additional 
amount of unbound stars that can escape detection in the 
present data. 

In particular, we use the approximations developed in 
D'Onghia et al. (2010) to compute the fraction of ICL that 
a disc galaxy with a dynamical mass equal to 5 x 10 11 Mq 
(of which 2 x 10 10 Mq is in the disc) produces passing on a 
parabolic orbit within the potential of a typical group. Such 
a galaxy has a stellar mass of the order of 10 10 Mq . Figure 5 
shows that the smaller the pericentric distance the larger the 
fraction of light stripped into the intergalactic medium. The 



18 In addition, we note that comparison of models and observa- 
tions is not straightforward even when considering simple stellar 
population models (Conroy & Gunn 2010 and references therein). 
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Figure 5. The fraction of ICL released by the model galaxy (see 
text) orbiting within the potential of a typical group computed 
as a function of pericentric distance (from D'Onghia et al. 2010). 
Due to the quasi-resonance phenomenon of the tidal interaction, 
the prograde encounters are more efficient in stripping material 
than the retrograde ones. 



fraction of ICL computed for the model galaxy is given by 
the ratio of the number of unbound stars to the total amount 
of stars bound to all the galaxies of the group. This galaxy, 
passing at 30-40 kpc from the group center, can contribute 
up to 3% of the total light in the group galaxies. Hence, 
slow tidal interactions likely provide a marginal contribu- 
tion to the total stellar budget of a low-mass fossil group 
as CXGG 095951+0212.6, but a non-neglible one for a fos- 
sil group as massive as CXGG 095951+0140.8. Given these 
rough estimates, we will estimate the ICL contribution to 
the total stellar luminosity of the two fossil groups under 
study in the future, thanks to the deep B- and R-band pho- 
tometry recently obtained by us with the Wide Field Imager 
(Baade et al. 1999), mounted at the Cassegrain focus of the 
MPG/ESO 2.2 m telescope in La Silla, Chile. 



5 CONCLUSIONS 

For the first time, the galaxy stellar mass function (GSMF) 
was computed, down to a stellar mass M stars ~ 10 9 Mq and 
in a robust statistical way, for two fossil groups discovered at 
2 « 0.4: CXGG 095951+0212.6 and CXGG 095951+0140.8. 
These groups have also well-determined total masses (M200) 

0' 



equal to 1.9 (±0.41) x 10" and 9.5 (±0.42) x 10" M, 
respectively. This was possible because of the high-quality 
multi-wavelength information available for the 2 deg 2 area 
of the Cosmic Evolution Survey (COSMOS; Scoville et al. 
2007) where the two systems were initially detected in X- 
rays (Finoguenov et al. in preparation; G09). In particular, 
the identification as fossil and the selection of the member- 
candidate galaxies of either group rest on robust photomet- 
ric redshifts, and iab ^ 22.5 mag limited spectroscopy. 
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Interestingly, the large-scale structure (LSS) in 
which the two fossil groups are embedded is very dif- 
ferent: CXGG 095951+0212.6 appears isolated, whereas 
CXGG 095951+0140.8 belongs to a LSS that covers the en- 
tire 2 deg 2 area of the COSMOS (corresponding to a cross 
size of about 25.5 Mpc at z « 0.4) and is traced by a total 
of 28 X-ray emitting groups. 

Comparison of the GSMF of these fossil groups and the 
composite GSMF of the COSMOS X-ray emitting groups at 
0.3 ^ z ^ 0.5 that span a similar range in total mass (Gio- 
dini et al. in preparation) highlights a lack of star-forming 
galaxies with M stars > 10 10 M Q in the two fossil groups. 
Consistently, the distribution of their photometric member 
galaxies in the b — r vs i colour-magnitude diagram reveals 
a well-defined red sequence and a lack of significantly bluer, 
luminous galaxies out to the virial radius -R200 of either 
fossil group. In addition, the total stellar mass fraction of 
these groups within 0.7i?2oo is relatively large for their to- 
tal masses (G09), in particular for the less massive one (i.e., 
CXGG 095951+0212.6). 

This evidence suggests that the overall conversion of 
(cold) gas into stars was accelerated and/or more efficient 
in fossil groups with respect to the average X-ray emitting 
group in the same mass range and at similar redshift. The 
star-formation activity in member galaxies rapidly ceased 
afterwards. At the same time, no significant infall of galaxies 
with M stars > 10 10 M Q took place in the last 3-6 Gyr for 
either fossil group at z « 0.4. 

If the processes shaping the galaxy component of a fos- 
sil group seem to be independent of the LSS, the progenitor 
of an observed fossil group may depend on it. On the basis 
of the measured stellar mass fractions (G09) and X-ray-to- 
optical luminosity ratios (Giodini et al. in preparation), we 
propose that compact groups are the most likely progen- 
itors of observed low mass and isolated fossil groups, like 
CXGG 095951+0212.6. Conversely, the relatively early infall 
of massive satellites likely originates massive fossil groups in 
dense regions of the LSS, like CXGG 095951+0140.8. 
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